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PACS 74.20.De - Phenomenological theories (two-fluid, Ginzburg-Landau, etc.) 

Abstract. - We propose a new stage of confiment-decofinment transition, which can be observed 
in laboratory. In two-gap superconductors (SCs), two kinds of vortex exist and each of them 
carries a continuously variable fraction of the unit flux quanta $o = hc/2e. The confined state of 
these two is a usual vortex and stable in the low temperature region of the system under a certain 
magnetic field above H c \ . We see an analogy to quarks in a charged pion. An entropy gain causes 
two fractional vortices to be deconfined above a certain temperature. 
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Introduction. — Topological defects furnish fascinat- 
ing problems in physics, and keep attracting a lot of con- 
cerns. In general, they lead the quantization of physical 
quantities. Vortex is one of the most important example of 
topological defects [1]. It is well known that in supercon- 
ductors (SCs), the requirement that the order parameter 
is a single- valued function guarantees the quantization of 
the magnetic flux of a vortex in the unit of <I>o = hc/2e. 

The fractionalization of topologically quantized num- 
bers have also been researched extensively in many dif- 
ferent contexts [2]. Vortices with fractional fluxes have 
been discussed in p-wave supcrfiuid 3 Hc [3] and e?-wavc 
SCs with a grain boundary [4] (See also a notice [5]). The 
fractionalization comes from unconventional structures of 
order parameters for these superfluid and SCs without any 
contradiction to their single-valued nature. 

Here we discuss vortex state in two-gap SCs [6] which 
has two superconducting gaps on two separate Fermi sur- 
faces. The two-gap superconductivity can be seen in a 
lot of SCs in solid state and has been paid special atten- 
tion. MgB2 is a recent typical textbook [7]. The possibil- 
ity of the two-gap superconducting state is also pointed 
out in 2H-NbSc 2 [8], PrOs 4 Sbi 2 [9], Y 2 C 3 [10], Ca-doped 
YBCO [11], and liquid hydrogen [12]. The fractionaliza- 
tion of the flux quanta has also been discussed in two-gap 
SCs [13,14]. Tanaka considered a wire ring and showed 
that a trapped flux in the ring can be an arbitrary fraction 
of $o [13]. Babaev considered a planer geometry and dis- 
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cussed the spontaneous fractional vortices associated with 
the Berezinskii-Kosterlitz-Thouless (BKT) transition [15] 
in the weak coupling limit of an internal Josephson inter- 
action with respect to the relative phase of the two gap 
functions [14]. 

Besides these pioneering works, searching for the ap- 
pearance of the fractional vortices under other setups is 
challenging and important since it gives new viewpoints 
of the problem and also new possibility for experimental 
discovery. In this paper, we propose a novel mechanism 
for the appearance of fractional vortices in the two-gap 
SCs. We consider the system under an applied magnetic 
field and show that a vortex with the unit flux is divided 
into two vortices with continuously variable fraction of <f>o 
by entropic effect. This phenomena can be seen as de- 
confinement of two fractional vortices. We see that in 
the confined state, two vortices behave like quarks in the 
charged pion [15]. Hence the present theory gives rise to a 
new stage of confmemcnt-decofinement transition, which 
can be examined in laboratory. 

Fractional vortices in a two-gap SC. Let us dis- 
cuss the property of fractional vortices in a type-II two-gap 
SC. We consider two gap functions A L and As (|A L | > 
| As |) opened on two different Fermi surfaces. Their pair- 
ing symmetry does not matter. We can consider not only 
s-wave but also unconventional pairing states [16]. We 
apply a magnetic field along z-direction to a sample of 
thin film geometry, i.e., d < mi7i(£g, £f ), where d is the 
sample thickness and and £f are the coherence lengths 
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of the gap functions along z-direction. The gap functions 
are then expressed as, Al,s(i") = | AL,s( r )| exp[— i6>L,s(r)], 
where r is the two dimensional vector in the xj/-planc. 
The Ginzburg-Landau (GL) free energy for two-gap sys- 
tem is straight extension of conventional one and given 
in Rcf. [17]. We consider the GL free energy in the Lon- 
don limit, namely, we neglect the spacial dependence of 
|Al(i")| and |As(r)| and obtain 
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The characteristic points of this free energy are the 
presence of two typical length scales, Xl,s(T) = 

(4e 2 /-r L ,s|A L ,s(T)|7?i 2 c 2 )~ 1/2 , where # L ,S are the coef- 
ficients of the gradient terms in GL free energy, and of 
the Joscphson-type interaction with a temperature depen- 
dent coupling T(T) = 27|A L (T)||A S (T)|. This coupling 
causes the locking effect for the relative phase 9^ — 9s- 
The equation of motion for the gauge field A is obtained 
from 5Fql/5A = as, 

A(T) 2 V x H = - A + ^V9 L + ^Vfl s , (2) 
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where 
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M(T) A 2 (T) 
is the London penetration depth, and 

A 2 (T) 
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The following boundary conditions arc imposed, 



V x V0 L , s (r) 



27rz<5 2 (r-R L ,s) 



$o and its ratio is determined by the ratio \l(T)/\s(T) 
(See, Eqs. (3), (4), and (5)). When A L (T)/A S (T) is 
limited to unity (such a limitation could come from the 
crystal symmetry, if we consider unconventional states in 
two dimensional irreducible representations in a system 
with single Fermi surface [16]), we have only the half 
flux $ L (T) = ®s(T) = $ /2, but there is no such a 
limitation in the present system and fractional vortices 
besides the half flux vortex are possible to exist when 
Al(T) 7^ Xs(T). It should be also emphasized that these 
fractional fluxes vary continuously with |Al(T)|/|As(T)|, 
which determines the ratio Xi,(T)/Xs(T). 

There is a possibility to observe these fractional vortices 
when their separation is large enough. It would be small 
in the low temperarture region due to the Josephson type 
interaction in the GL free energy (1), but is possible to 
be large in the high temperature regime in order to gain 
entropy. To check this scenario, we calculate the energy 
cost and the entropy gain for a certain configuration of 
two vortices with $l(T) and <f>s(T). 

Configuration energy. — We can eliminate the 
gauge field from Eq. (1) by using the equation of motion 
(2), and obtain the energy cost, 
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{V(0 L - S )Y -T(T) cos(0 L - 6 S )] , 

where K{T) = {$ L (T)$ S (T)} /16tt 3 A(T) 2 . The vortex 
solution in the thin film geometry for the single gap SC 
has been investigated [18]. The extension to the two-gap 
case is straightforward. We introduce the effective London 
penetration depth in the thin film A cff (T) = \ 2 (T)/d [18]. 
When A c ff (T) is comparable to the system size, we obtain 



" 6i{T)-2nK(T)ln 11 



i=L,S 



Aeff(T) 



£rel,(9) 



(6) where 



which would be relevant for the dilute vortex system with 
H ~ H c \. The singlevaluedness of two gap functions is 
guaranteed except for Rl and Rs because of the 27r wind- 
ing. From Eq. (2), we obtain the London equation for the 
two- gap case 



A 2 (T)V x V x H + H 
$ L (T)z5 2 (r-R L ) + $ s (T)z<5 2 (r 



Rs 



(7) 



This equation tells us important implications. When 
Rl = Rs = R, the equation reduces to the usual Lon- 
don equation A 2 (T)V x V x H + H = $ z<5 2 (r - R), 
and we have a vortex with the unit flux passing through 
R. When Rl 7^ Rs, however, we have two vortices with 
$ L (T) through R L and <J>s(T) through R s . It should 
be emphasized that <&l(T) and $s(T) are fractions of 



e L , s (T) = {$ 2 ! s(r)/167r 2 A c ff(T)}ln(A c ff(T)/a,s(T)) 

are the creation energy for L- and S- vortices, ^l,s{T) = 
hvp' S /|Al,s(7 1 )| the coherence lengths of L- and S-gap 
functions, each of which should be in the same order of 

L s 

each vortex core, and v F ' the Fermi velocities. The log- 
arithmic term in the r.h.s is the magnetic interaction be- 
tween two vortices, which is repulsive and shows the log- 
arithmic behavior. The last term has a sine-Gordon form 
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When 7 = 0, r(T) = and £ r ci oc In |R L - R s | [14]. Such 
a case can be applied to the superconductivity in liquid 
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Fig. 1: A draw with four excitations of the string. The tick 
line denotes the string connecting two fractional vortices with 
a distance R. 



hydrogen [12]. We go beyond the weak coupling limit of 
7, which would be significant for the two-gap supercon- 
ductivities in solid state, e.g., MgB2 [7]. 

In Eq. (10), the spatial dependence of the relative phase 
looses the elastic energy. The Josephson-type term favors 
to be 8 L - 9 S = 2nn for 7 > and 8 L - 8 S = (2n + 1)7T for 
7 < 0. But we have the boundary condition (6). There 
should be a line- like singularity between separated vortices 
associated with a 2-7T jump of the relative phase. This 
"string" is the two-dimensional extension of the kink in 
the sine-Gordon model in one dimension [19]. Therefore, 
for large 7, 

E rc \ = e st (T)L st + const., (11) 

where e st (T) = d^J K{T)\T{T)\/tt + 0(1/7) is the string 
creation energy per a unit length and L st > |Rl — Rs| is 
the length of the string. To minimize H v , L st — > |Rl— Rs|- 
This indicates that there is a linear confinement poten- 
tial between L- and S-vortices. We have also the two- 
dimensional Coulomb type repulsion. Then, two fractional 
vortices is analogous to massive " quarks" with " fractional 
electric charge" $l and <&s, except for the spatial dimen- 
sionality [15]. The L-S confined state corresponds to the 
"charged pion" (total charge +e of pion corresponds to 
the unit flux $0 = + &s of L-S confined state). 

Entropy, confinement, and deconflnement. 

Let us examine the thermodynamic stability of the sep- 
arated state with a distance R = |Rl — Rs| = N£ ^ 0, 
where £ = max(^L, £s) an d we refer it as the cut-off length 
of this system at the short length scale. We introduce an 
artificial square lattice with a constant £ so that two vor- 
tices sit on one of the lattice axes. We discuss the free en- 
ergy F = H V — TS, where H v is given in Eqs. (9) and (11). 
There are two contributions to the entropy S = S\ oc + 5 s t , 
the former one comes from the location of two vortex cores 
Sioc — 2fceln(ri/£ 2 ) (ft; the area of the system) and the 
latter the string configuration S s t- 

The lowest energy configuration of the string is the 
straight line connecting two vortices, where L st = R. The 
kink-like configurations of the string would be thermally 



Fig. 2: A dependence of tds 
deconflnement temperature. 



Tde C /T c , where Tdec is the 



excited and each excitations contributes £ to L st (See, Fig. 
1). As shown in Fig. 1, there are two kinds of excitations 
one of which is directed to the "upward" and the other 
" downward" . The numbers of the upward and downward 
excitations are equal each other so that the string is ended 
by two fixed points (vortices) and we denote them as n. 

We consider states with L s t ^ R, i.e., n <C N. 
By using the dilute approximation we obtain S s t = 
k B In N C n N-nC n ~ k B (SR/0 [1 - hidR/2R], where SR = 
2n£. We find the value of SR which minimize F at fixed R 
is SR (R) = 2i?exp [— £,e s t/k-BT]. Let us define a function 
F(R) = F(R,SR {R)). We introduce a{T) as the value 
of the vortex distance R at the mimimum of F{R) and it 
shows a diameter of the confined region of two fractional 
vortices. We obtain 



a(T) 



2TrK(T)/e st (T) 



2k B T 



5(T)e at (T) eX P 
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Remind that the solution of 1 —2/x- exp[— x] = is x = 
LambertW(2) ~ 0.85. Then, a(T) -> 00 [20] when 



i(rKtcn 

LambertW(2) 



0. 



(13) 



The temperature that satisfies Eq. (13) would correspond 
to the deconfinement temperature and we refer it as Td ec - 
We emphasize that it is surely to be Td ec < T c , where 
T c is the gap opening temperature. The reason is that 
£(T)e st (T) monotonically decreases to increase T and van- 
ishes at T c . Then, the l.h.s. of Eq. (13) changes its sign 
once between the zero temperature and T c . This state- 
ment is crucial since vortices cannot exist above T c and 
the argument becomes inconsistent when Td ec > T c . 

We can also show the presence of Td ec when we con- 
sider states with L st >> R. Let us remind that the num- 
ber of the configuration N con = n n for the n-step self- 
avoiding random walk, where \i is a constant and equal 
to 3 for the 2D square lattice case. For the string with 
sufficiently large L st , N con would be well approximated 
by the random walk result and we obtain N con ~ /i 
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Table 1: Parameters which would be relevant for MgB2 [7] 



f(0)/A(0) 


0.32 


7 


1.6/cV/cell 


d 


10 A 


X 


0.33 


a 


2.2 


k 


2.1 



Then, the entropy of the string is proportional to L st [15] 
and the string free energy F s t — T s t(T)L s t, where the co- 
efficient T st (T) = {e B t(T)-kBT/i(T)-\n/j,} denotes the 
string tension and is positive in the low temperature re- 
gion, i.e., fractional vortices are confined. Above a certain 
tcmpctarure, T s t(T) < 0, namely, the string looses its ten- 
sion and dcconfinemcnt occurs. Td ec is obtained from the 
equation T st {Td ec ) = 0, which is equivalent to Eq. (13) 
when we take /j, = exp [LambertW(2)] ~ 2.35. 

To obtain Tdec, we should know the temperature de- 
pendence of two gaps. By using GL equation for two-gap 
superconductor with sufficiently large 7 [17], we check that 
|A L ,s(T)|/|A L , s (0)| = /c % /T^ r t, where t = T/T c and k is 
a constant which can be obtained numerically. Then, Eq. 
(13) becomes A\J\ — t/t = LambertW(2), where, in the 
case of Ks/Kl = 1 [17], 
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4tt 5 / 2 \ A(0) J V (x + 1/x) 



and x = I A s (0)|/| A L (0)| . We see from Fig. 2 that t dec = 
Tdec/Tc depends on A strongly when A ~ 0(1). By using 
the parameters in Table 1 which would be relevant for 
MgB 2 , we obtain A ~ 5.7 and t dec ~ 0.97. 

We discuss BKT transition temperature associated with 
the usual vortex (6*l = #s) obtained from /cbTbkt = 
$o/327r 2 A e //(T B KT) (See, Eq. (1)). By using the param- 
eters in Table 1 , we see £bkt — 0.90, which is lower than 
tdec we have obtained. The present argument for tdec is 
relevant for the dilute vortex system but does not include 
the multi-vortex effect which would be important above 
Tbkt, since a lot of thermally activated vortices and anti- 
vortices exist above Tbkt even in the system under the 
applied magnetic field considered here [21]. Then, the re- 
sult would be corrected in this case. It should be noted 
that the coefficient A includes £(0) explicitly, on the other 
hand, the equation for Tbkt does not. So, Tdec can be 
lower than Tbkt in the system with sufficiently small £ (0) 
and the present argument can be precise. It is well known 
that high-T c cuprates have rather small coherence length 
in the order of 10 A. Therefore a two-gap superconductor 
Ca-doped YBCO [11] seems to be one of the best candi- 
dates to apply the present argument. 

We comment on the relation to the renowned decoupling 
transition [22] . The renormalization group (RG) approach 
of the 2D sine-Gordon model discussed in Ref. [22] sug- 
gests that the cosine term in Eq.(10) becomes irrelevant 
above a certain temperature T*, which indicates that the 
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Fig. 3: a(T) for A = 1.6. In this case, t dec ~ 0.81. Note that 
the unit of the vertical line is £(T) = £(0)/\/l — t. Here, we 
use the relation a(0) ~ 0.3£(0) obtained from Table 1. 

decoupling between #l and 0g occurs. The lower bound of 
T, is obtained from T* = 8irK(T*)d in the limit 7^0 [22] . 
This RG analysis takes into account string loops (closed 
strings) only, i.e., it is assumed that the relative phase 
#l — #s is single valued. The open string (See, Fig. 1) 
should be taken into account since it would be created by 
the thermal fluctuations. Actually, we have shown that 
the open string is stabilized by the entropic effect and 
the dcconfinemcnt occurs at Td ec - For sufficiently small 
A, Tdec can be lower than the lower bound value of T», 
i.e., the deconfinement can occur even in the temperature 
region that the cosine term is relevant. 
We also find that 



a(T) 
a(0) 



1 



2t 



Ay/T^t 



exp {-A/1 - t/t} 



(14) 



where a(0) = 27rK - (0)/e 8t (0). The plot of a(T) is given in 
Fig. 3 for A — 1.6. We see that even in the confined phase 
it seems to be possible to observe each of fractional vor- 
tices in the region of large a(T) by using some appropriate 
vortex imaging techniques. 



The authors are grateful to J. Akimitsu, E. Babaev, 
G. Baskaran, G. Blatter, H. Fukuyama, M. Hayashi, I. 
Hcrbut, T, Ishida, R. Ikeda, K. Ohashi, and Y. Tanaka for 
their critical comments. This work is financially supported 
by Grant-in- Aid for Scientific Research from Japan Society 
for the Promotion of Science under Grant Nos. 15540370, 
16740226, and 18540381. 

REFERENCES 

[1] See, for example, P. G. de Gennes, Superconductivity 
in Metals and Alloys (Perseus Books, MA) 1999. 

[2] See, for reviews, A. J. Niemi and G. W. Semenoff, 
Phys. Rep., 135 (1986) 99; R. B. Laughlin, Rev. Mod. 
Phys., 71 (1999) 863; G. E. Volovik, The Universe in a 
Helium Droplet (Oxford University Press, Oxford) 2003 



T..A 



Dcconfinement of Vortices with Continuously Variable Fractions of the Unit Flux Quanta in Two-Gap Superconductors 



[3] G. E. Volovik and V. P. Mineev, Pis'ma Zh. Exp. Teor. 
Fiz., 24 (1976) 605 [JETP Lett, 24 (1976) 561] 

[4] M. Sigrist, D. B. Bailey, and R. B. Laughlin, Phys. 
Rev. Lett, 74 (1995) 3249; D. B. Bailey, M. Sigrist, 
and R. B. Laughlin, Phys. Rev. B, 55 (1997) 15239 

[5] The fractional vortex is defined as a vortex with a flux less 
than $o /2 in Ref. [4] , while less than $o in this paper. 

[6] H. Suhl, B. T. Mattias, and L. R. Walker, Phys. Rev. 
Lett., 12 (1959) 552 

[7] T. Muranaka et al., Frontiers in Superconducting Mate- 
rials, edited by A. V. Narlikar (Springer- Verlag, Berlin) 
2005 

[8] T. Yokoya et al, Science, 294 (2001) 2518 

[9] G. Seyfarth et al, Phys. Rev. Lett., 95 (2005) 107004; 

M. YOGI et al, J. Phys. Soc. Jpn., 75 (2006) 124702. 
[10] S. Akutagawa et al, (unpublished). 
[11] J. H. Ngai et al, cond-mat/0703423. 
[12] N. W. Ashcroft, J. Phys., A129 (2000) 12. 
[13] Y. Tanaka, J. Phys. Soc. Jpn., 70 (2001) 2844; Phys. 

Rev. Lett, 88 (2002) 017002. 
[14] E. Babaev, Phys. Rev. Lett, 89 (2002) 067001; Nucl. 

Phys. B, 686 (2004) 397; E. Babaev, Asle Sudbo and 

N. W. Ashcroft, Nature, 431 (2004) 666. 
[15] J. B. KOGUT, Rev. Mod. Phys., 51 (1979) 659. 
[16] M. Sigrist and K. Ueda, Rev. Mod. Phys., 63 (1991) 

239 

[17] M. E. Zhitomirsky and V.-H. Dao, Phys. Rev. B, 69 

(2004) 054508 
[18] J. PEARL, App. Phys. Lett, 5 (1964) 65 
[19] T. H. R. Skyrme, Proc. Roy. Soc, A247 (1958) 26. 
[20] Actually, a(T) would have a cut off determined by the 

inverse of the square root of the applied magnetic field. 
[21] S. Doniach and B. A. Huberman, Phys. Rev. Lett 42, 

1979 (1169) . 

[22] J. P. Rodriguez, J. Phys. Cond. Mat, 9 (1997) 5117; 
Phys. Rev. B, 62 (2000) 9117; R. Goldin and B. 
Horovitz, ibid., 72 (2005) 024518; S. Coleman, ibid., 
11 (1975) 2088; S. E. Korshunov, Europhys. Lett, 11 
(1990) 757. 



